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We address the issue of identifying the mediators of effective interactions in cuprates supercon¬ 
ductors. Specifically, we use inversion theory to analyze Raman spectra of optimally and over-doped 
La 2 - :c Sr a ;Cu 04 samples. This allows us to extract the so-called glue function without making any 
a priori assumption based on any specific model. We use instead two different techniques, namely 
the singular value decomposition and a multi-rectangle decomposition. With both techniques we 
find consistent results showing that: i) two distinct excitations are responsible for the glue function, 
which have completely different doping dependence. One excitation becomes weak above optimal 
doping, where on the contrary the other keeps (or even slightly increases) its strength; ii) there is a 
marked temperature dependence on the weight and spectral distribution of these excitations, which 
therefore must have a somewhat critical character. It is quite natural to identify and characterize 
these two distinct excitations as damped antiferromagnetic spin waves and damped charge density 
waves, respectively. This sets the stage for a scenario in which superconductivity is concomitant 
and competing with a charge ordering instability. 

PACS numbers: 74.72.-h, 71.45.Lr, 74.20.Mn, 78.30.-j 


I. INTRODUCTION 

The issue of the mediators of high-temperature super¬ 
conductivity in cuprates, the so-called “glue issue”, is still 
unsolved. On the one hand, there is the early proposal of 
an essentially instantaneous magnetically mediated pair¬ 
ing forming incoherent singlets eventually condensing in 
the superconducting state.® On the other hand, it has 
been proposed that slow, nearly critical collective excita¬ 
tions mediate a retarded pairing. In this latter framework 
overdamped spin waves are a natural candidate due to 
the proximity of the superconducti ng phase to the Mott 
insulating antiferromagnetic phase PS How strong is the 
instantaneous or the retarded character of the magnetic 
supercondu cting glue is likely just a matter of quantita¬ 
tive balance!^ On the contrary, there are other propos¬ 
als suggesting that the superconducting mediator could 
also have a non-magnetic character. In this regard, pri¬ 
mary candidates are charge-ordering fluctuations, which 
were proposed as a source of retarded glue already long 
agoEHU This last proposal has steadily acquired impor¬ 
tance with the progressive increasing experimental evi¬ 
dence of charge ordering as a common feature of cuprates 
via neutron scatte ring] ^ NMR and NQRp® EXAFS,^ 
ARPESp-®E] stm]^®P and Raman scattering J2IH23] un¬ 
til it has become a compelling evidence due to the recent 
direct observation by resonant X-ray scatteringP3f22] 

The evidence of charge ordering and the simultane¬ 
ous natural presence of spin waves emanating from the 
nearby AFM state raises the issue of coexistence, in¬ 
terplay, and possible causal relationship between charge 


and spin degrees of freedom. An early point of view 
is that charge ordering arises from Coulomb-frustrated 
charge segregation due to short-range (essentially instan¬ 
taneous) magnetic inter action d 31 * 32 ' or/and due to non¬ 
magnetic (e.g., phononic^ 3 2H, charge-transfer®!,...) at¬ 
tractive forces. In both cases, the concomitant presence 
of charge-ordering fluctuations “enslaves” the spin degrees 
of freedom, allowing them to survive up to large doping. 
This tight entangling of spin and charge degrees of free¬ 
dom is also an intrinsic feature of the ‘stripe” concept.®^ 
A complementary, more recent, view is that collective 
retarded spin fluctuations may give rise to attraction in 
the partic le-hol e channel inducing a charge-density wave 
instability.^®^ At the moment, it is unclear and debated 
whether these two points of view are really distinct or 
just some kind of “egg-and-chicken” issue arising from 
an underlying continuous interplay between charge and 
spin degrees of freedom. In any case it seems by now 
rather natural to assume that charge-density waves and 
spin-density waves coexist and it is quite interesting to 
investigate the relative weight of these different excita¬ 
tions in mediating the effective electron-electron interac¬ 
tion in cuprates in the different doping and temperature 
regimes. This is precisely the focus of the present work. 

The presence of two distinc t sou rces of glue in 
cuprates was found both in opticaP^i, ARPE^^l anc i 
Rauiaipd-DJ experiments. These experiments identify the 
so-called “glue function” a 2 F(co), introduced by Eliash- 
berg to characterize the spectral distribution and the 
strength of the electron-phonon coupling in the supercon¬ 
ducting pairing.®! From optical conductivity it is found 
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that the glue function has a double structure with a 
rather narrow and temperature dependent peak at low 
frequencies (below 10 3 cm -1 ) and a broad structure ex¬ 
tending up 3000 — 4000 cm -1 . From the characteristic 
frequency ranges it is quite natural to attribute the first 
peak to phononic excitations, while the mediators re¬ 
sponsible for the broad peak should be diffusive mag¬ 
netic excitations.^ However, neutron scattering experi¬ 
ments show that the m agnet ic excitations have a rather 
rich spectral structure which, besides the broad 
high-frequency peak, also display a low-frequency peak. 
Moreover, time-resolved optical spectroscopy finds that 
also the low-frequency peak could be related to elec¬ 
tronic excitations Pil Thus the attribution and identifi¬ 
cation of the various excitations responsible for the ef¬ 
fective electron-electron interaction is still an open rel¬ 
evant issue. In this framework, the analysis of Raman 
spectra is interesting because this technique is able to 
explore selectively different regions of the electronic Bril- 
louin zone allowing to extract more detailed information 
on the electronic excitations than the optical conductiv¬ 
ity. This selectivity of Raman spectra was already ex¬ 
ploited in Ref. 1231 to identify and characterize the double 
nature of the electronic mediators, but the analysis was 
carried out fitting the experiments by analytic expres¬ 
sions arising from a specific assumption on the form of 
the charge and spin collective excitations. In particular 
the standard diffusive form of nearly critical overdamped 
collective modes near gaussian quantum critical points 
was assumed. This not only allowed a good fitting of the 
spectra, but allowed to attribute the markedly different 
shape of the Raman spectra in the different channel (the 
Big and B 2 g scattering channels obtained by differently 
polarizing the incoming and outgoing photons of the Ra¬ 
man scattering) to different, spin and charge, excitations, 
which acted differently on the spectra according to their 
different characteristic wavevectors. Although success¬ 
ful, this approach was starting from specific assumptions 
on the form of the mediators, which had the single spec¬ 
tral structure of overdamped modes. It is therefore quite 
important to test and deepen the conclusions of Ref.[231 
with a more flexible and general approach. In this work 
we apply two different numerical inversion techniques to 
extract the glue function from the Raman spectra and we 
draw consistent conclusions: i) two distinct excitations 
are responsible for the glue function, which have com¬ 
pletely different doping dependence; ii) there is a marked 
temperature dependence of the weight and spectral dis¬ 
tribution of these excitations, which therefore must have 
a somewhat critical character. Indeed, if they were just 
the result of low-energy phonon-mediated interaction and 
high-energy local excitations due to strong-correlation ef¬ 
fects, any temperature dependence would hardly occur. 

Our analysis is purposely carried out on optimally and 
overdoped samples to highlight the specific features of the 
bosonic interaction mediators only. This analysis could 
of course be extended to the underdoped regime, where, 
however, the presence of a pseudogap in the fermionic 


spectrum introduces additional temperature, doping, and 
momentum dependencies, which would be superimposed 
to those of the bosonic mediators. This would make the 
interpretation of the data much less conclusive and infor¬ 
mative. 

The paper is organized as follows: In Section II we 
describe the two numerical methods adopted to extract 
the glue function from the Raman spectra. Sections III 
and IV contain the results of our analysis of the glue 
functions and of their low-frequency behavior. Finally, 
Section V reports and discusses our concluding remarks. 
Appendices A and B contain technical details on the im¬ 
plementation of the inversion procedure, while Appendix 
C contains details on the extraction of the low-frequency 
properties of the glue function. 


II. GLUE FUNCTION AND RAMAN SPECTRA 


In this work we analyze Raman spectra from Ref.[23l 
in the Bi g and B 2g channels in La 2 _ a; Sr a; Cu 04 (LSCO) 
samples for different temperatures and doping levels. 
Electronic Raman scattering is a bulk (nearly surface- 
insensitive) probe and it measures a response function 
x(z) analogous to that of the optical conductivity. The 
Raman response function can be expressed^ as a func¬ 
tion of the so-called memory function M(z) as x( z ) = 
XoM{z)/[z + M{z)\ for complex arguments z, where Xo 
is the (real) “bare” response function computed in ab¬ 
sence of any scattering process. The imaginary part of 
the response function, x"■> can t> e written as: 

X"H = uM"{u) 

Xo + w] 2 + [M"(u ])] 2 ' U 

Here, the real and imaginary parts of the memory func¬ 
tions, M(u>) = are function of the real 

frequency w, and they are related by a Kramers-Kronig 
(KK) transformation. In turn, M"(u>) can be expressed^ 
in terms of the glue function a 2 F(z) by means of the in¬ 
tegral expression 



dz K(uj,z)a 2 F(z), 


with the kernel 
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Therefore, x" is a non linear functional of a 2 F. The 
extraction of a 2 F from experimental data is possible by 
approximating it as a linear combination of suitable ba¬ 
sis functions as a 2 F{oj) « c a^a(w). Once a ba¬ 

sis is chosen, one may optimize the coefficients c a so as 
to fit the experimental data. In this way, Eq. ([2]) can 
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be rewritten as M"(u) = ]r a CgA"(w), with ^a( w ) = 
dz K(u, z)(j> a (z), so that Eq. 0 becomes 

X"(u,{c<*}) = _ EgD^gM _ 

XO Eg C a A' a (uj) + w] 2 + Eg Cgv4"(w)] 2 ' 

(4) 

where A' a (uj) is the KK transform of A"(w). The co¬ 
efficients c Q are then used as fit parameters. We used 
and compare the results obtained with two very differ¬ 
ent choices for the basis functions, the first one based 
on a multi-rectangle decomposition (MRD) of the glue 
function, and the second based on the singular vector 
expansion of the kernel K. 

• Multi-Rectangle Decomposition (MRD). The 

glue function is approximated by a piecewise con¬ 
stant function, which corresponds to choosing a 
partition of the frequency axis oji < ... < w Q < 
... < uin+i and N non-overlapping box functions 
<t> a {u) = 1 for u} a < lo < w „+i and 4> a (u) = 0 
elsewhere, as the basis functions. The fitting pa¬ 
rameters c a are the heights of the box functions. 
In the MRD approach the functions A" can then 
be computed analytically (see Appendix [A] for fur¬ 
ther details). 

• Singular Vector Decomposition (SVD). The 

starting point of this approach is the expansion of 
the integral kernel Eq. (J2]) as a sum of diadic oper¬ 
ators: 

OO 

K(w, z) = ^2 cr a 4> a (u)(fi a {z) 

0=1 

where the singular values a a are nonnegative and 
in decreasing order, and {ip a } and {<f> a } are sets of 
orthogonal functions, called the left and right sin¬ 
gular vectors, respectively. In this case we approx¬ 
imate the glue function using as basis the first N 
right-singular vector s of th e singular vector decom¬ 
position of the kerneP^kMI. This allows us to define 
A"(w) = cr a ip a (oj). Notice that the integral kernel 
of our physical problem has to be properly regular¬ 
ized before the expansion thus leading to slightly 
different expressions; we refer to the Appendix [B] 
for a detailed derivation. 


In both cases coefficients c a are computed by mini¬ 
mizing the square distance A 2 between the theoretical 
response function and the experimental one. In order to 
do this, the functions A! and A" are computed at the 
experimental points u) t , i = 1,..., iV exp . We then look 
for 


JV, : 


min A 2 ({c q }) = min ^ [x"xp,i - Xth,i({cg})] , (5) 

{c “ } {c “ } i= i 


where Xth ?:({ c «}) — X"(Ai- { c g}) is computed according 
to Eq. Q, and x" xp>i = XexpE) are the measured val¬ 
ues. The minimization is constrained by the requirement 


a 2 F(z) > 0, which can be written in the MRD case as 
c Q > 0, and in SVD case as J2 a c a (^a{z) > 0, for all pos¬ 
sible values of z. Further details about the choice of the 
functions <j> a and xo> and in the discretization and in the 
fitting procedures can be found in Appendices [0 (MRD) 
and [B] (SVD). 


III. B lg AND B 2g GLUE FUNCTIONS: DOPING 
AND TEMPERATURE DEPENDENCE 


We analyze Raman spectra of La 2 _ a; Sr a ,Cu 04 sam¬ 
ples focusing on three different doping levels (x = 
0.15, 0.20, 0.25) and three different temperatures (ap¬ 
proximately 50, 100 and 200 K; the true temperatures 
for different doping levels are slightly different), to high¬ 
light the doping and temperature evolution. In Fig.[T] 
we show the glue functions obtained with both MRD 
and SVD techniques by analyzing data of three differ¬ 
ent doped samples (x = 0.15, 0.20, 0.25) at T ~ 50 K, 
for both the Bi g and B 2 g channels. 


The glue function obtained by SVD has an oscillatory 
character, due to the fact that it is obtained by sum¬ 
ming up oscillating functions with increasing frequency, 
similarly to a Fourier expansion. (See Fig. B1 for repre¬ 
sentative examples of the basis functions.) On the other 
hand, the MRD approach may yield very irregular glue 
function. Some sort of indeterminacy is an inevitable 
consequence of the structure of the problem. Integral 
kernels act as low-pass filters, and therefore only the low- 
frequency components of the glue function, related to the 
largest singular values, can be reliably determined from 
the (noisy) data. (See Fig. |B3| for the magnitude of the 
singular values of the kernel.) Instead of relying on arbi¬ 
trary smoothing techniques, we proceeded in a pragmatic 
twofold way. On the one hand we considered integrated 
quantities, which are relatively insensitive to both discon¬ 
tinuities and oscillations in the glue functions, and there¬ 
fore allow to extract robust relevant informations. On 
the other hand, one can notice that despite the different 
“looks” of the glue functions obtained with the two ap¬ 
proaches, some features are common and therefore more 
reliable. Indeed, both approaches yield in both channels 
two structures: a low-frequency peak (w < 500 cm” 1 ) 
and a broader structure in the range 500-3000 cm” 1 . The 
Big channel appears to be suppressed as the doping is 
increased, while the secondary structure in the I? 2 g glue 
function shifts to lower frequencies as the doping is in¬ 
creased. These behaviors are also present at T = 100 
and 200 K. This latter case is reported in Fig. [2] Due 
to the behavior of the kernel at high temperature (see 
discussion in Appendix B 4), the structures of the glue 
functions are rougher, but still the main features found 
at T = 50 K persist, namely a generic two-peak structure 
and a different doping dependence of the two channels. 

We can quantify these features by computing the total 
spectral weight Wtot of the glue function and the fraction 
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Figure 1. (Color online)Fitted glue functions a 2 F(z) using the SVD (red) and MRD (black) methods, for three different doping 
levels [(a) and (d), x = 0.15; (b) and (d), x = 0.20; (c) and (f), x = 0.25], for both B\ g [(a)-(c)] and B 2 g [(d)—(f)] channels at 
low temperature (~ 50 K). The insets show the experimental data (grey circles) and the response functions computed from the 
fitted glue functions. 


of the weight below a cutoff frequency w = 500 cm 1 , as 

POO 

Wtot = / a 2 F(z) dz , 

W Q _ £ a 2 F(z) dz 
Wtot / 0 °° a 2 F(z) dz 

In particular, the ratio between W& and W tot gives us the 
magnitude of the low-frequency structure with respect to 
the higher-frequency one. Results obtained by MRD and 
SVD are shown in Fig.[3| At all doping level and all tem¬ 
peratures we find good quantitative agreement between 
the integrated quantities computed via MRD and SVD. 
Data show little dependence on temperature, while the 
effect of doping is much stronger. As described before, 
the total area of the R lg channel decreases markedly as 
the doping increases, with the total spectral weight at 
x = 0.25 being half the weight at x = 0.15. On the other 
hand, the fraction of weight below 500 cm -1 appears to 
be roughly constant, suggesting that the two structures 
are proportionally suppressed. The B 2 g channel glue 
function shows quite the opposite behavior. In fact, the 
total weight is only slightly dependent on temperature, 
while the fraction of weight below uj = 500 cm -1 increases 
from ~ 10% at x = 0.15 to ~ 20%. 


IV. LOW-FREQUENCY ANALYSIS: 

TEMPERATURE AND DOPING DEPENDENCE 

The previous analysis show little temperature- 
dependence of the global properties of the glue function. 
However, the low-frequency part of the Raman response 
function has a clear dependence on temperature. A low- 
frequency expansion of the glue function allows us to 
study the behavior of the glue function at frequencies 
w < T. As detailed in Appendix |C 1[ the relation be¬ 
tween the Raman response function and the memory 
function M can be approximated for finite temperature 
and small frequencies as 


x"M = Xo 


UJ 

K 


+ 0(u 3 ) 


( 6 ) 


Therefore, the low-frequency behavior of the response 
function is governed by the imaginary part of the memory 
function at zero frequency, Mg. Interestingly, Mg' can 
be directly connected to the low-frequency behavi or of 
the glue function a 2 F(z). As detailed in Appendix |C 2[ 
Mq is proportional to an “effective” slope s e g of the glue 
function, obtained by a weighted average of a 2 F(z)/z 
with a fast decaying weight, similar to exp {—z/T). The 
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Figure 2. (Color online) Fitted glue functions a 2 F{z) using the SVD (red) and MRD (black) methods, for three different doping 
levels [(a) and (d), x = 0.15; (b) and (d), x = 0.20; (c) and (f), x = 0.25)], for both B\ g [(a)-(c)] and B 2g [(d)—(f)] channels at 
high temperature (~ 200 K). The insets show the experimental data (grey circles) and the response functions computed from 
the fitted glue functions. 


relationship between s e ff and Mq is given by 

3Mq 

SeS ~ 4t r 3 T 2 ' 


(7) 


If the glue function is approximately linear for z < T, 
then s e ff is approximately equal to its slope. On the other 
hand, if the glue function has two peaks, one for z < T 
and the other one in z > T, the quantity s e s is related to 
the magnitude and steepness of the first peak only, sim¬ 
ilarly to the ratio Wq/W tot discussed before. As already 
discussed, while the derivative of the glue function can¬ 
not be reliably estimated from the fitted glued function 
without some kind of smoothening procedure, integrated 
quantities, like s e ff in ([7]), are much better suited for the 
analysis of the glue function properties. 

We find that Mq displays a clear temperature trend, 
which is well described as a T 1 / 2 scaling for both the 
channels. Furthermore, the two channels have a marked 
opposite dependence on the doping. If we use a simple 
linear model Mq oc (a + bx) to include the effect of the 
doping, we obtain the following expressions: 


Mq(T, x, B lg ) « (190 - 570 x) ■ VT cm^K' 1 / 2 , 
Mq(T, x, B 2g ) « (300 x) ■ Vf cm -1 K -1/2 . (8) 


These scalings reflect into similar scalings for the slope of 
the Raman response at zero frequency and effective slopes 
s e ff of the glue functions (see Appendix [C] and figures 
therein for more details). These results are in agreement 
with the analysis shown in the previous section. The 
total weight W to t in the B\ g channel decreases almost 
linearly with the doping, but not W^/Wtot, therefore 
reducing s e ff. On the other hand, in the B 2g channel, the 
total weight is constant, but HAoo/I'I'tot increases with x; 
as the low-frequency peak “drains” spectral weight from 
the second one, the slope s e ff also increases. 


V. DISCUSSION AND CONCLUSIONS 

The Raman spectra in the two, Bi 9 and B 2 9 , channels 
are markedly different: the B\ g spectra have a rounded 
rise at low frequency until a new rounded growth takes 
place above 2000 cm” 1 , while the B 2g spectra have a flat- 
tish shape at low/intermediate frequencies until again a 
rounded growth takes place above 2000 cm -1 . This dif¬ 
ference naturally stems from different physical mecha¬ 
nisms ruling the low/intermediate-frequency scattering 
mechanisms. According to the most popular scenarios 
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Figure 3. (Color online) Panels (a) and (b): Total area 
(weight) Wtot = f dz a 2 F{z) under the glue functions as 
a function of doping x. Panels (c) and (d): weight frac¬ 
tion below 500 cm -1 , computed as Wmo/Wtat with W 500 = 
fg°° dz a 2 F(z). SVD results are shown in panels (a) and (c), 
while MRD results are shown in panels (b) and (d). Shades of 
red: Bi g channel. Shades of blue: Bi g channel. Temperature 
increases from dark to light colors. 


of superconducting cuprates, spin and charge collective 
excitations are the most likely candidates as mediators 
of scattering. In a previous work 23 we noticed that the 
different direction and magnitude of the characteristic 
wavevectors of these excitations produce different con¬ 
tributions in the two channels of the Raman spectra. 
Specifically, it was found that charge fluctuations with 
characteristic wavevectors q ( , sa (±7 t/ 2, 0), (0, ±7r/2) (in 
reciprocal lattice units) are more visible in the low- 
frequency spectra in the £? 2 g channel. On the con¬ 
trary, spin fluctuations with characteristic wavevectors 
q s ~ (±7r =F x, ±7r), (±7t,±7t =f x) (in reciprocal lattice 
units) contribute more in the low-frequency spectra of 
the Big channel. At higher frequencies both modes con¬ 
tribute to both channels. While the previous analysis as¬ 
sumed a phenomenological quantum critical form of the 
spin and charge excitations leading to analytic but ob¬ 
viously model-dependent glue functions, our main goal 
here is to put under scrutiny the previous results with an 
unbiased general analysis. Therefore within a memory- 
function approach we extracted the glue function directly 
form the Raman spectra without any assumption. To this 
purpose, we used two inversion methods: the SVD and 
the MRD. The overall shapes of the extracted glue func¬ 
tions are similar with both methods (see Fig.JlJ: in both 
channels they show a rather narrow low-frequency peak 
below 500-1000 cm -1 and a broad hump up to 3000- 
4000 cm -1 . Despite the similarity of the two channels, 


a closer inspection shows that the doping dependency is 
just opposite, with the B\ g glue function rapidly decreas¬ 
ing with doping above the optimal, while the glue func¬ 
tion extracted from the Big spectra stays roughly con¬ 
stant. More than this, one can notice that the frequency 
above which the Bi g glue functions vanishes shrinks pass¬ 
ing from the 3000-4000 cm -1 range, at x = 0.15, to about 
2000 cm -1 , when x = 0.25. This points to the natu¬ 
ral interpretation that spin waves underly the B\ g glue 
function at all frequencies, while charge excitations are 
responsible for the low/intermediate-frequency scatter¬ 
ing in the i? 2 g channel. In this channel spin waves also 
contribute, but mostly at high frequencies. When dop¬ 
ing suppresses the spin-mediated scattering, the overall 
Big glue decreases as well as the high-frequency part of 
the B 2 g . The marked difference of the doping depen¬ 
dence in the two channels is also visible from both the 
low-frequency and total spectral weights of the glue func¬ 
tions (see Fig. [ 3 ]). A marked increase of the low-energy 
excitations involved in the Bi g channel is also clear form 
Fig. 0r.dk 

The glue functions extracted here are in overall agree¬ 
ment with those obtained from other techniques. A 
peaked feature below 1000 cm -1 is also present in op¬ 
tical experiment 39 -— and i n low-energy spin e xcitat ions 
revealed by inelastic neutron scattering (INS) 47 ^!). At 
higher energies (up to 0.3-0.4 eV) broad humps are 



Figure 4. (Color online) Main panel: Temperature depen¬ 
dence of the inverse of the effective slope of the glue function 
at low frequency, l/s e ff, for different channels (Bi g , circles; 
Bo g , triangles) and doping levels ( x = 0.15, yellow; x = 0.20, 
pale green; x = 0.25, dark green). Flere, s e ff has been com¬ 
puted from the Mq fitted using the SVD method. Similar re¬ 
sults are obtained with the MRD. Inset: The same data in a 
log-log scale, highlighting a seemingly T 3 / 2 behavior (dashed 
line). 
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found in optical spectra like those obtained within our 
analysis, which can be identified as due to spin exci¬ 
tations both from neutron scatterm;®! and from reso¬ 
nant inelastic X-ray scattering (RIXS)^^. At the same 
time the re is in creasing evidence from both NMRpU and 
RIX^SHzHUSni that CDW are ubiquitous low-energy ex¬ 
citations in cuprates. It is therefore natural to identify 
them with the excitations involved in the B 2 g channel. 

Our results also display a general agreement with 
theoreticaP-B and experimental^, ARPES results in 
B^S^CaCi^Og+a samples. In particular there is full 
consistency concerning the coexistence of two distinct 
collective modes and with the softening of the mode(s) 
near the “hot” antinodal region. It is noticeable that 
this agreement is present despite the obvious differ¬ 
ences between the Raman response and the electron self¬ 
energy. In particular the former is resolved with respect 
to the electron momentum, while it is integrated over all 
the momenta of the bosons dressing the single-particle 
fermionic propagator. On the other hand, the Raman 
response obeys some (partial) cancellations when the bo¬ 
son momenta and the Raman form factors interplay 2 ^. 
From Ref. 323 one sees that the mode, which softens in 
the antinodal region, gets more strongly coupled in the 
underdoped region and could be identified with our spin- 
fluctuation mode. The (momentum averaged) energy of 
the other mode stays more constant over the electronic 
momentum space and it acquires relatively more weight 
increasing the doping above optimal doping. This leads 
to the identification of this mode with our mixed phonon- 
CDW mode, which becomes prominent when doping in¬ 
creases above the optimal level. Concerning the tem¬ 
perature dependence, INS find below T c a sharpening of 
the low-energy part of the spin-fluctuation spectra (see, 
e.o.l43j). accounting for the clear observation of this mode 
in the ARPES spectra (Ref.s 01] and FT21) below T c (the 
contribution of the other broad featureless, nearly tem¬ 
perature independent part of spin fluctuations being nat¬ 
urally buried in the broad incoherent parts of the elec¬ 
tron self-energies). The softening of the charge modes 
extracted (mainly) from the Raman (i? 2 g) spectra is only 
expected to occur near the antinodal region and therefore 
in ARPES may be overshadowed by the softening of the 
other spin mode. Our selectivity in momentum allows 
instead to disentangle the contributions of the modes to 
the different Raman channels and to detect a substantial 
softening of the charge mode as well. 

Our work focuses on temperatures above T c and finds 
a distinct temperature behavior for the overall and the 
low-energy parts of the spectra. Specifically, we notice 
a rather weak dependence of the overall spectra, which 
becomes more pronounced in the low-frequency weight. 
This indicates that the excitations might have a low- 
frequency component with a marked temperature depen¬ 
dence. Indeed, as explained in Sec. IV, exploiting an¬ 
alytic low-frequency expansions of the memory function 
and of the glue functions extracted, e.g., via the SVD 
method, we are able to identify a substantial tempera¬ 


ture dependence of the the imaginary part of the memory 
function at zero frequency. This is related [see Eq. ©1 
to the slope of the Raman response near uj = 0. In turn, 
this is related to the effective slope of the glue function 
s e ff defined in Eq. |7]), which is a rather robust quantity 
extracted from a weighted average of a 2 F{z)/z that is 
only weakly sensitive on the oscillatory character of the 
SVD fitted glue functions. These strong temperature de¬ 
pendencies at low frequency are suggestive of some form 
of critical behavior of the mediating excitations. If for a 
while we borrow the analytic expressions of the glue func¬ 
tions calculated in Ref. 123[ we find m/g ~ l/s e ff, where m 
is the mass of a collective mode (i.e., the minimum energy 
needed to excite it) and is proportional to the square of 
the inverse correlation length of the corresponding fluctu¬ 
ations, while g represents the coupling between fermionic 
quasiparticles and the collective modes. Thus, although 
in this work we purposely avoided any assumption on the 
mediator, we interestingly find that our results support a 
marked temperature dependence of its mass, as it is seen 
in Fig. [4] 

Within a quantum critical scenario, the masses are 
expect to vary linearly with T in the quantum critical 
regime, and saturate in the quantum disordered regime. 
Of course drawing any conclusion on this dependence 
with three temperatures only is out of question, but a 
clear indication of nearly critical mediators, with a cer¬ 
tain tendency of the mass to flatten as a function of tem¬ 
perature at low T, can still be obtained. 
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Appendix A: MULTI-RECTANGLE 
DECOMPOSITION 

In this case the strategy to minimize the functional 
in Eq. ([5]) is based on a histogram decomposition of the 
glue function, the fitting parameters being the heights 
of the histograms, in the N different frequency inter¬ 
vals (bins) in which the frequency axis is partitioned, 
and the prefactor \o (that also include all multiplicative 
parameters needed to match the theoretically calculated 
response with the measured one). The fitting procedure 
searches for a local minimum, beginning from a starting 
guess, and u sing the direct search algorithm of Hooke 
and Jeevesf^Sir 

Using as starting point Eq. ([2|, we divide the integra- 






tion range into N non-overlapping intervals as follows: 


N 


'M = ^2 / dzK(u},z)a 2 F(z), 


(Al) 


OL= 1 ' 


where j = 1,... ,N, and to avoid confusion with the fre¬ 
quencies u>i at which the Raman response is experimen¬ 
tally measured, we called here a a and b a the frequencies 
identifying the a-th bin, with a a < a Q +i = b a . In every 
integration interval of Eq. (Al I, the glue function is to be 
considered a constant (the height of the bin) and can be 
taken out the integral and considered as a multiplicative 
factor indicated with c Q , so we can rewrite the equation 
ill) in a matrix form 


N 


M>,) = E c “ AM ", 


(A2) 


with 


AM" a E f dzK(u)i,z). 

J a a 


The integra tion domain in Eq. © extends between 0 and 
oo. In Eq. (All we have taken b jy = 8000cm -1 = z max 
assuming that the integrand vanishes at higher frequen¬ 
cies. The minimum frequency (acting as a natural cut¬ 
off within our procedure) was typically taken as ai = 
10 cm -1 = z m i n . Rather than adopting a homogeneous 
mesh, with constant b a — a a , we adopted a logarithmic 
mesh, with constant b a /a a , to enhance the sensitivity at 
low frequency. The value of the integral in the bin with 
extremes [a, 6] appearing in Eq. (Al), with a, b > 0, is 


A M"(u) = 2TO~ 1 a 2 F(a) ( 2Q In 

l smh(a) 

+ — l j ) — ^{b -\- cj ) zi(b — oj ) 

1 

— —2 [A(2 ka + 2kuj) — A(2 ka — 2 kO)\ 

k=1 

V . 

+ —2 [A(2 kb + 2kui) — A(2 kb — 2kui)] |, 

k =1 


where a = a/(2T), b = b/(2T), a; = w/(2T), 5(x) = 
x[ln(2sinh|x|) — ||a:|], A(x) = e - l x lsign(a’), and the last 
two lines contain rapidly convergent series. For numerical 
reasons, we approximated ln(sinh |x|) rs |a;| — ln(2) when 
|ar| > 40. 

Once the imaginary part of the memory function is ob¬ 
tained, the real part is achieved through the KK trans¬ 
formation 


1 r + °° M"(z) 

M\w) = -P / clz-bi. 

7T / ^ Z-tJ 


(A3) 


Since the M" (uj) funct ion is known at discrete points, to 
calculate the integral (A3) the integrand function is re¬ 


placed with a continuous broken line obtained by joining 


the heights of the bin with straight lines. So the analyt¬ 
ical form of M"(u>) to be included in (A31 will look like 
M"(z) = Az + B, in which the coefficients A and B are 
calculated within each bin. To calculate the integral (A31 
we have subdivided the integration interval in subinter¬ 
vals which are delimited by the same point where M"(oj) 
is known. Apparently, we can have two cases: u> is or is 
not an extreme of integration. In the first case the results 
of the integral (A3) for the subinterval delimited by, say, 
[ui \, W 2 ] is straightforwardly 


A(uj 2 — kq) T (Auj -(- R) In 


u >2 — u> 

UJ\ — OJ 


while in the latter case two neighboring bins [wi,w — £] 
and [w + ^, W 2 ], treating the divergence in the sense of the 
principal value, have to be considered. When the results 
for the two neighboring bins are summed, the divergent 
part cancels and the same expression is obtained. 

We have to remark that, in general, multiple choices 
of a 2 F yield the same memory function M", due to 
the presence of nonzero solutions to the matrix equation 
A M” a c a = 0. It is therefore impossible to fit an 
unique glue function a 2 F without imposing additional 
constraints (see, e.g., Ref.G>nj). We restrict the possible 
glue functions by imposing the following constraints: 

1. a 2 F(u;) > 0, and therefore c a > 0; 

2. The initial guess for the glue function in the min¬ 
imization procedure has vanishing spectral weight 
for uj > 4000 cm -1 , assuming that the collective 
modes live at lower frequencies. 

Such constraints do not completely avoid the uncertainty 
about the glue function. However, the relevant features 
of the calculated glue function are robust enough, upon 
varying the doping and the temperature, as discussed in 
Sec-Hill 

We performed our calculations comparing the results 
obtained with N =8,12, and 25 bins, to ensure that the 
main features of the glue function were robust with re¬ 
spect to variations in the MRD scheme. 


Appendix B: SINGULAR VECTOR 
DECOMPOSITION 

In this section we describe in a detailed way the SVD- 
based fitting procedure. The first step is to discretize 
the kernel K(u>, z) in Eq. ([3|. This is done by using two 
(possibly different) meshes {wi} and {zj}. Logarithmic 
meshes, or combinations of linear (at small frequencies) 
and logarithmic (at large frequencies) meshes, can be 
used. The integration measure is associated to the kernel, 
so that Eq. ([ 2 ]) becomes: 


M" = ^2 K ijOc 2 Fj . (Bl) 

3 
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The SVD, although computationally demanding, pro¬ 
vides the best low-rank approximation of the kernel K, 
a result known as the Eckart-Young theorem. In this 
sense, it is the ideal tool to approximate solutions of in¬ 
tegral equations. However, a naive application of the 
SVD to the kernel in Eq. (Bl) produces basis vectors cj> a 
with undesirable properties. (1) On physical grounds, 
we require that the glue function goes to zero at large 
frequency (at least as \/z, see Ref.lK2l). However, the in¬ 
tegration kernel K is long-ranged, in the sense that it is 
nonzero for large u and 2 . As a result, the singular vec¬ 
tors (f>a(z) do not go to zero for large z, making hard to 
constrain a 2 F to be small at large frequencies. (2) Even 
worse, the kernel K(co,z) diverges as ~ 1/z for 2 —> 0; 
therefore, the SVD of the kernel produces basis vectors 
(j> a which also divergent in 0. This is in sharp constrast 
with the fact that the glue function should go to zero 
at small frequencies in order to ensure the convergence 
of the integral, Eq. ©• As described in Appendix [Aj 
within the MRD approach these two problems are eas¬ 
ily solved by choosing bins with lower edges larger than 
zero (thus providing a cutoff to the integral) and upper 
edges smaller than a maximum frequency ~ 8000 cm' 1 . 
In the SVD context, the solution to these problems re¬ 
quires to expand a properly regularized kernel instead of 
the original one; we will detail our procedures in the next 
sections. 


1. Large frequencies 


In order to perform numerically the integral in Eq. ([2]) , 
some cutoff frequency fl has to be introduced. Due to the 
long-range nature of the kernel, the <f> eigenfunctions are 
strongly dependent on the cutoff. Therefore, the choice 
of the cutoff should be motivated by physical reasons, 
and we have to check that the results are only weakly 
dependent on the particular cutoff choice. 

Let us note that the integration kernel K(uj,z ) goes 
to zero for z > uj + O (e~ ( - z ~ ul ^ T ) [see Eq. |3|]. This 
means that high frequency (z > fl) components of the 
glue function do not contribute to the memory functions 
at low frequency (uj < fl) and, vice versa, that only the 
low-frequency part of the glue function can be fitted from 
low-frequency data. 

We can explicitly set to zero the high-frequency part 
of the glue function introducing a cutoff function q(z) as 

a 2 F(z ) = q(z)f(z), (B2) 

where q(z) is 1 at low frequencies and goes to zero at high 
frequencies, and / does not diverge in the z —> oo limit. 
Many choices for the cutoff function q(z) are possible, for 
instance a power law decay 


«(*) 


1 


or an exponential decay, e.g., 

«<*> = \ 

for suitable parameters fl, 7 and A. We find that the 
results depend only weakly on the cutoff, as the sup¬ 
port of the fitted glue functions is concentrated before 
2000-3000 cm -1 . Therefore, we will write M"(ui) = 
f dz K(uj, z)q(z)f(z) and fit f(z) to the experimental 
data using the cut off kernel K(u>, z)q(z). After a func¬ 
tion f(z) has been determined, the corresponding glue 
function is obtained from Eq. (|B2[). 



2. Small frequencies 

The divergence in 2 = 0 of the kernel implies that 
the right eigenfunctions (which are used as a basis func¬ 
tion for the glue function) are divergent in zero. Since 
a 2 F( 0 ) = 0 , we want to approximate the glue function 
using a set of functions for which <fi a ( 0) = 0. This can be 
accomplished by expanding a suitable regularized kernel 
instead of the original one. We can write: 

K(uj, z) = K r (uj, z) + K d (z) , (B4) 

where K R is regular in z —> 0, and I\d(z) = 47 rT )z. We 
can use the SVD decomposition for the regular part of 
the kernel, and study the effect of the second term on the 
solution. Of course, Kd(z) is still divergent, but the right 
eigenfunctions of K R go to zero fast enough to regularize 
the integral. 


3. Rewriting the integral equation 


Putting together the cutoff definition, Eq. (B2), and 
the decomposition into regular and divergent part, 
Eq. (B4), we are able to recast the initial equation Eq. ([ 2 ]) 
in the following form: 

M"(w) = fdz [.K R (uj,z)q(z)\f(z ) 

+ j dz[K D (z)q(z)\f(z). (B5) 

We then perform the SVD on the new kernel K’(uj, z) = 
[K r (uj, z)q(z)\, rather than on the original kernel 
K(u,z), as follows: 


N 


K r (uj, z)q(z) « Y v a il’a(u)<l> a (z), 


(B 6 ) 


where N is the number of singular vectors we use to ap¬ 
proximate the kernel. We can hence rewrite Eq. (B5l by 
expanding f(z) = J2 a c a (j) a (z), obtaining 

M"(z) = Y c «K(z), 


1 + (z/flp ’ 


(B3) 
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Figure Bl. (Color online) Singular vectors from the SVD ex¬ 
pansion of the regularized kernel, Eq. (B6l, at T = 100 K 
for some representative values of a = 1,2,4, 8. (a) ip a (id) 

functions (b) Kramers-Kronig transform ip a {oj) of the 
functions (c) (p a (z) functions. The eigenfunctions are quali¬ 
tatively similar for all temperatures considered in this work. 


where 


A"(w) = a a tp a (u}) + Sa a , 

= J dz K D (z)q(z)<t> a (z). 

The real part of the memory function, M' (z) can be ob¬ 
tained as a function of the coefficients c a by taking the 
KK transform of the functions A"(w) (or, equivalently, 
of the ip a (in) functions). Eq. is then obtained by plug¬ 
ging the expressions for M' and M" into Eq. ©• 


4. Fitting procedure 


A simultaneous fit of both yo and the expansion co¬ 
efficients {c a } is often unstable, with Xo being pushed 
towards very large values. Therefore, we decided to fit 
the coefficients c a at fixed yo, and to systematically study 
the effect of different xo- As explained in the main text, 
the coefficients {c a } which fix the glue functions are ob¬ 
tained by the minimization of a function A 2 ({co,}), i.e., 
the square distance A 2 between the theoretical response 
function and the experimental one. 

It is useful to introduce the primitive of the glue func¬ 
tion W(z) = fg dz'a 2 F(z'), and W max = IV(z max ), with 
Zmax being the largest frequency in the z-mesh. There¬ 
fore, the quantity 1 — W {z)/W (z max ) equals the fraction 
of the area under a 2 F at frequencies larger than z. 

\ shows an example (data corresponds to the B lg 


B2 


Fig. 

channel, T = 105K, doping x = 0.25). As one can see 
in panels (a) and (b), small values of yo do not allow for 
a good fit, while for yo > 17 the A 2 is generally small. 
However, we see that using a large yo produces glue func¬ 
tions with a sizable spectral weight at high frequencies, 



Figure B2. (Color online) (a) Fit of x" using different yo 
values. We show the goodness of fit, as measured by A 2 
(see Eq. 0 ). and the ratio between the area under the glue 
function for lj > 2500 cm -1 to the total glue function weight 
Wmax- Best fit is obtained for yo ~ 17. (b) Smaller values 
of yo do not yield a good fit, with y " h being systematically 
smaller than the experimental susceptibility, (c) The fitted 
glue function increases its total area, in particular at high 
frequencies, as yo is increased. The lines shown in panel (b) 
[or (c)[ correspond to the best fit value yo = 17 (blue thick 
line) and the yo values shown in panel (a), smaller (or larger) 
than 17 (thin lines). All points have been computed using the 
same mesh (N w = 536, N z = 507) and cutoff (power-law, see 
Eq. ( |B3| ) with Q = 3000 cm -1 and 7 = 4 ). 


as it is clear from panels (a) and (c). In this case the 
choice yo = 17 gives the best combination of fit quality 
and reduced spectral weight at high frequencies. 

We finally address the question of how to choose the 
number of singular vectors. The contribution of each 
singular vector <j> a to the susceptibility is mediated by 
the singular values cr a . The number of singular vectors 
which have to be kept into account depends on how fast 
a a decreases with the index a. As one can see in Fig. |B3| 
the decay of the singular vectors strongly depends on 
temperature. At high temperature, the kernel can be ap¬ 
proximated by just using a few terms of the SVD Hence, 
the number of singular vectors used to approximate the 
kernel should decrease with T, too. Moreover, we find 
that the convergence of the fit gets considerably worse if 
too many singular vectors are used. A reliable choice is 
using 40 SVs for T ~ 50 K, 25 SVs for T ~ 100 K, and 15 
SVs for T ~ 200 K. 


Appendix C: RELATING LOW-FREQUENCY 
BEHAVIOR OF RAMAN RESPONSES AND 
GLUE FUNCTIONS 


1. Raman response and memory function 


Some informations on the glue function can be directly 
extracted from the low-frequency behavior of the exper¬ 
imental Raman response function y". In fact, let us ex¬ 
pand both the real and imaginary parts of the memory 
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Figure B3. (Color online) (a) Magnitude of the singular values 
a a of the kernel /\r(oi, z)g(z), at three different temperatures, 
normalized by the value of the first SV. (b) Number of singular 
values a a larger than 10 -6 x ai as a function of temperature. 
In both panels the same mesh (N w = 536, N z = 507) and 
cutoff (power-law, see Eq. (B3l with SI = 3000 cm -1 and 7 = 
4) was used. 


2.5 1 - 


m 1-5 


X Using M" 0 from glue functions 

Q Using fitted M' 0 (T,x) 


0.5 ^ 
0.5 


ztX 


O x X 


✓ o< 
x o ^ 


1.5 

Modeled x" at 50 cm 1 


c 


2.5 


function near zero as 


M'(u) = M[oj + 0(w 3 ), 

M"{u) = M[! + M 2“ 2 + 0(w 4 ), 


where we exploited the opposite parity of the functions 
M' and M". Substituting into Eq. 0 we obtain 


Xo 


Ml' 


MS 


[(M[ + l) 2 - M"M"] +0(w 5 ). 

(Cl) 


Therefore, the low-frequency behavior of x is controlled 
by the imaginary part of the memory function, computed 
at zero frequency. As described in the main text, Mg' is 
found to have a marked dependence on temperature and 
doping, as described by Eq. ([8]) and shown in Fig. Cl 



Figure Cl. (Color online) The red points are MS using the 
SVD procedure, as a function of temperature and doping, for 
the two channels. The surfaces represent the fits Ml'(T,x) = 
VT’(a + bx). The surface parameters are a = 190 and b = 
—570 for the Bi g channel, while they are a = 0 and b = 300 
for the B 2 g channel, where both a and b are in K - 1 / 2 cm -1 
units. 


Figure C2. (Color online) Comparison between the experi¬ 
mental value of the Raman response functions at u = 50 cm -1 
and the values obtained using the approximation x”( w ) ~ 
Xotu/Mo', with Mg' being directly computed from the (SVD-) 
fitted glue functions of with the fitted models in Eq. | 8 }. In 
both cases, xo is set to the values obtained during the fitting 
procedure. 


The expansion Eq. (Cl) allows an estimate of the Ra¬ 
man response function x /, ( w ) i n th e low-frequency limit. 
Fig. |C2| shows how the experimental x"(o;) computed at 
u) = 50 cm -1 (the lowest nonzero frequency available for 
all the datasets) compares with the first order expansion 
Xow/Mq, where Mg' can be either the value computed 
by integrating the fitted glue functions, or the value ex¬ 
tracted from the fits in Eq. (|8|. The agreement is far 
from perfect, although the correct trend is respected. We 
have to remark that the single experimental points used 
in the comparison can be strongly influenced both by 
the limited number of data points at low temperature, 
by random instrumental errors and/or by the presence 
of phonon excitations. As a consequence, the scaling re¬ 
sults, Eq. ([8]), obtained using Mg (arising from global fits 
to the data) are much more consistent than the ones ob¬ 
tained with x" alone. 


2. Memory function and glue function 


Equation ([2 1 can be used to relate Mg' and the glue 
function as follows 


MS = 


d z K(0, z)a 2 F(z ) 


= 2tt 


dz —cosech 2 ( —) a 2 F(z) . 
2 T V2 TJ 


The kernel A'(0, z) has a 1/z divergence at z = 0 and dies 
exponentially as exp(— z/T). Therefore, the support of 
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o) (cm ') z (cm ’) 


Figure C3. (Color online) Low frequency behavior of Raman 
response function and glue functions a 2 F at T ~ 50 K. 
The Big channel is shown in panels (a) and (b), while the 
£> 2 g channel is shown in panels (c) and (d). The straight 
lines in panels (a) and (c) correspond to x" ~ Xo^/Mg, see 
Eq.§; the lines in panels (b) and (d) correspond to a 2 F(z) = 
s e ffZ, where s e ff is the effective slope computed from Mg , see 
Eq.0. All the lines are computed with the Mg and the xo 
extracted using the SVD fitting procedure. 


K(0,z) is essentially the interval 2 : £ [0, X 1 ], so that Mg' 
is a measure of the spectral weight of the glue function 
at frequencies u> < T. In particular, at low temperature, 
the slope at zero frequency of the glue function is the 
dominant contribution in the integral. This suggests the 
definition of the quantity 


•Seff — 


3M" 
4tt 3 T 2 ‘ 


(C2) 


It is easy to show that, in the T —> 0 limit, s e ff coincides 
with the slope da 2 F(z)/dz at z = 0. This can be seen 
by expanding in series the glue function around z = 0 

as a 2 F(z) = C\Z + c-iz 2 + C 3 Z 3 + _ The coefficients 

Cj depend in principle on both temperature, doping and 
channel. Then, we can compute the contribution of each 
term to Mg' (this may not be a convergent series, it de¬ 
pends on how fast the coefficients c n decay with n; how¬ 
ever, we assume that it is at least an asymptotic series 


for low enough T): 

pOO 

<-£/ K(0,z)z n dz 

n> l Jo 

pOO 

= S^2n(2T) n+1 / cosech 2 (j/) y n+1 dy . 

n> 1 

The function J(x) = J 0 °° cosech 2 (t) t x dt is well defined 
for complex-valued arguments in the domain Re(x) > 1, 
and can be expressed as J(x ) = 2 1_a: £(x)r(x + 1), 
Re(x) > 1, where T(x) and £(x) are Euler’s gamma and 
Riemann’s zeta functions, respectively. The first few val¬ 
ues of the integral are J(2) = 7r 2 /6, J(3) = 3£(3)/2 and 
J(4) = 7 t 4 /30. This leads to the result 


M" 


^2 tt (2 T) n+1 J(n+1) 


n> 1 


4t tT c n T n x C [n + l)T(n + 2) 


n> 1 

~ %T 2 cx + 24C(3)T 3 c 2 + ^T 4 c 3 + ... (C3) 
3 15 

Therefore, s e g ~ C\ if the first term dominates the se¬ 
ries in Eq. (C31. Using the known value of the integral 


J(2), we see that another way to define s e ff, equivalent 
to Eq. (C2|, is 


^eff 


r d- 


2 F(z) 


(it ) 2 cosech 2 (2^) 


fo°° d “(#) 2 cosech 2 (^) 


This expression represents an average of the 
quantity a 2 F(z)/z with a fast-decaying wei ght 
« {z/T) 2 exp(— z/T) for z T. We show in Fig. C3 
how the values Mg' computed from the SVD-fitted glue 
functions can be used to estimate the low frequency 
behavior of the Raman response function x"(w) [panels 
(a) and (c)[, and the effective slope s e ff. We see that 
the extrapolated slope of fits quite well with the 
experimental data, even when the data points are scarce. 
The effective slope of the glue function, instead, is 
systematically smaller than the slope of the fitted glue 
functions at z = 0. This is due to the glue functions 
being concave. However, s e g appears to be roughly pro¬ 
portional to da 2 F/dz\ z ~Q. In particular, this suggests 
that the dependence of s e ff and da 2 F’/d 2 :| z= o on doping 
are very similar. Fig. [4] in the main text reports the 
temperature dependence of the inverse effective slope 
s e ff, highlighting the strong temperature dependence of 
the low-frequency glue function, possibly related to the 
quantum critical behavior of the interaction mediator. 23 
Although the plot is limited to three temperatures, this 
quantity seems to be quite compatible with a power-law 


dependence s 


eff 


T 3 / 2 , as it is seen from the inset of 


Fig. [4] calling for further investigation. 
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